Establishing the Presence of Coherence in Atomic Fermi Superfluids: Spin-Flip and 
Spin-Preserving Bragg Scattering at Finite Temperatures 
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We show how the difference between the finite temperature T structure factors, called S-, associated with 
spin and density, can be used as a indication of superfluidity in ultracold Fermi gases. This observation can be 
exploited in two photon Bragg scattering experiments on gases which undergo BCS- Bose Einstein condensation 
crossover. Essential to our calculations is a proper incorporation of spin and particle number conservation laws 
which lead to compatibility at general T with two /-sum rules. Because it is applicable to general scattering 
lengths, a measurement of S- can be a useful, direct approach for establishing where superfluidity occurs. 
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One of the most intriguing aspects of the ultracold Fermi 
gases in the strongly interacting, attractive regime is the diffi- 
culty they pose for establishing superfluidity. The absence of 
indications from the density profiles, along with the presence 
of a normal state excitation gapare all responsible for this dif- 
ficulty. Most experiments fljj— 14£| which determine whether a 
given state is above or below the transition at T c , make use 
of a fast sweep of the magnetic field into the low field Bose 
Einstein condensation (BEC) limit, where the superfluid sig- 
natures are more straightforward. Thermodynamical proper- 
ties do not establish whether the phase is ordered or not, al- 
though they can determine where T c lies It would be 
particularly useful then to find a neutral fluid counterpart to 
the Meissner experiment for Fermi gases with attractive inter- 
actions. More specifically, we want to find a property which 
reflects the superfluid order parameter such that it is non-zero 
below T c and vanishes above. Of particular interest is Bragg 
spectroscopy [Q] which is capable of detecting true superfluid 
coherence. This is to be contrasted with radio frequency spec- 
troscopy [8] which is most suited for establishing pairing. 

It is the purpose of this paper to address this goal by pre- 
senting a theory of the spin preserving and spin flip dynam- 
ical responses at general temperatures in a homogenous gas. 
Both of these can, in principle, be measured via two photon 
Bragg spectroscopy. Our emphasis here is on the contrasting 
physics of spin density and particle density which we demon- 
strate allows a separation of coherent order from pairing ef- 
fects. While we focus on unitary gases, we also present re- 
sults associated with both the BCS and the BEC sides of reso- 
nance. Indeed, there has been considerable theoretical interest 
in combined studies of spin and particle density Bragg scatter- 
ing at zero temperature in BCS-BEC crossover theory II 9l l ldfl 
as well as of spin-flip Bragg experiments above T c II ill . 

Two aspects are essential to a proper theory of Bragg scat- 
tering, (i) At all temperatures the identities based on conser- 
vation laws for particle number and spin need to be satisfied: 
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where Sc,s(w, q) are the dynamical structure factors for the 
particle density (labelled C) and the spin (labelled S). Here 
Uf is the total density of particles. The first equation is the 
well known /-sum rule in the form associated with the gen- 
eralized BCS Hamiltonian. In past work El sum rules were 
found to be difficult to implement in the presence of pairing 
correlations, (ii) Below T c , collective mode physics in the 
charge channel is responsible for Eq. ©, not Pauli blocking 
as is sometimes claimed; this ensures current conservation and 
restores the U(l) gauge symmetry. 
We define 

S ± ((J,q) = ^[Sc(uJ,q)±S S (u;,q)] (3) 
From Eq. ([TJ an important additional sum rule then follows 

/oo 
du)U)S-(u,q)=0. (4) 
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The quantity 5_ (ui, q) represents a measure of "spin density 
and particle density separation". In the literature H the quan- 
tity is defined in the same way as S-. Importantly, S'-j-j,, 
has been argued 13, UM to be proportional to the density den- 
sity correlation function (p^(r)p±(r)). The sum rule of Eq.© 
indicates that Sfi times frequency integrates to zero. In strict 
BCS theory this sum rule is satisfied because 5_(w,q) (or 
equivalently Sfi(io, q)) has negative weight below T c and is 
identically zero above T c . In this paper we show that the be- 
havior found in strict BCS is general. On this basis, we argue 
that the association of the density- density correlation function 
(with opposite spin) and the quantity S^^(uj,q) (or equiva- 
lently S-(uj, q)J is very problematic. Collective mode effects 
are, in part, the reason that naive linear response theory (relat- 
ing (pt(r)pj,(r)> to S n ) fails below T c . 

The theoretical approach used here is based on the BCS- 
Leggett ground state wavefunction, which has been shown J^] 
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Figure 1 : The diagrams contributing to the full electromagnetic ver- 
tex A. Here the wiggly lines represent the T-matrix t pg and the 
dashed line represent the singular "condensate propagator" t sc . In 
AL2, A' is the full gauge invariant EM vertex. 



to give a reasonable description of the T = spin and den- 
sity dynamical structure factors as computed using quantum 
Monte Carlo simulations. Here we generalize this treatment 
of BCS-BEC crossover to finite T and must choose the appro- 
priate diagram set to respect the conservation laws in Eqs.([T|i 
and ||2). The present approach is most well suited to address 
these sensitive issues which arise in gauge invariant electro- 
dynamics; here in contrast to other crossover schemes IU3I1 
the superfiuid density is well behaved and monotonic in T. 
We build on the analysis outlined in Ref. IU4I1 which focused 
principally on collective mode effects. In this extension of the 
ground state to finite T, the pairing gap A(T) has two con- 
tributions associated with condensed (sc) and non-condensed 
(pg) pairs such that A 2 (T) = A 2 JT) + Ap g (T). The latter 
vanishes at T = and the former vanishes at T c and above. 



For spin-preserving and spin flip Bragg processes 11 1 511 the 
external field Hamiltonian contains contributions of the form 
/ <Pr[\ c J,^ a , + h.c] where \ G aa , oc <W and X s aa , cx 
8 a a' ■ Here ^(ifr) is the fermionic creation (annihilation) op- 
erator, a =t, I, t= — I and a = —a. While there would 
seem to be a similarity in these two channels we stress that 
there is also a physical difference. In the superfiuid phase, 
due to coherent singlet pairing, the response functions of spin 
and particle densities are different, whereas one may antici- 
pate that they are the same in the normal phase. Indeed, linear 
response theory based on the above Hamiltonian is incom- 
plete; it will not lead to a consistent treatment of the density 
preserving response, because it ignores collective mode ef- 
fects below T c . For the unperturbed Hamiltonian of the BCS 
type, and for singlet pairing, spin is conserved. However, be- 
cause the spin degrees of freedom are not associated with a 
spontaneous symmetry breaking, they do not lead to collec- 
tive phase modes upon condensation. 

In the density channel, linear response theory is more com- 
plex; we introduce a four-vector formalism with "electro- 
magnetic" (EM) field where \i = 0,1,2,3. The in- 
duced EM current J M is given by J"(Q) = K^"{Q)A V {Q), 
where, Q = = (w,q) is the four- momentum, and the 
EM response kernel, K^ v , can be written as K^ V (Q) = 
K% V {Q) + 8K» V {Q), where 8K^ V is the correction due to 
collective mode effects (3. Here 8K^ V , which is essen- 
tial for charge conservation (q^K^ v (Q) — 0), has a van- 
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Figure 2: (Color online) (a) Static structure factor (in units of np/2) 
at unitarity in the density channel S c (q) /2 (solid lines) as a function 
of q for two temperatures and S- (q) (dashed lines), (b) Temperature 
dependence of (q) showing that it behaves like an order parame- 
ter, as in BCS theory and for all Here kp, cf, and Tf are the 
Fermi momentum, Fermi energy, and Fermi temperature of a non- 
interacting Fermi gas with the same particle density. 



ishing denominator ifTill associated with the order parame- 
ter collective mode dispersion. In this letter, we are inter- 
ested in the 00-component of the EM kernel, which corre- 
sponds to the gauge invariant density-density response func- 
tion Xpp(^,q) = K 00 (u>, q). This, in turn, is related to 
the dynamic structure factor in the density channel, given 
bySc(w,q) = -icoth(^)Imx P p(o;,q), where x PP (w,q) 
represents the properly gauge invariant form of the density- 
density response function. In a similar way, the dynamic spin 
structure factor is S,s(u, q) = — icoth(^)Imxss(w, q), 
where xssi^, q) is the spin-spin response function. 

Now we turn to finite temperatures where there is little prior 
work addressing Bragg scattering in the crossover scenario. 
Charge conservation and gauge invariance which can be ex- 
plicitly shown to hold at the BCS level must be respected 11411 
even in the presence of pairing correlations. The best way to 
insure these is to establish that the /-sum rules Eq. (HJ and the 
related Eq. (O are satisfied 111 611 ■ For our BCS-Leggett ground 
state, extended to T ^ 0, the diagrammatic corrections to 
the EM vertex have been discussed [14] and are indicated in 
FiglD There are two types of contributions shown as Maki- 
Thompson (MT) and Aslamazov-Larkin (ALi 2 ) diagrams. 

The "bare" part of the EM response kernel is given by 
K v {Q)-^g^{l-g^) = 



2Y,^(K,K- Q)G(K)G(K - Q)A"(K, K-Q), (5) 

K 

where the bare and full EM vertex are X^(K,K + Q) — 
(1, i(k - |)) and A M = A p + SA SCfl + (5A pgA1 respec- 
tively. Here m is the particle mass and g^ v is the metric 
tensor. Note that the full vertex A does not include collec- 
tive mode physics and, thus, is not gauge invariant below T c , 
while the vertex A' contained in the AL2 diagram is gauge 
invariant. Importantly, consistency with the Ward identities 
will lead a cancellation B14I1 between the MT and AL terms 
|(ALi + AL 2 ) + MTfj; = 0. Throughout we use the super- 
script "C" to refer to the density response and the counterpart 
"S" for the spin counterpart. 

After imposing the Ward identities B1411 . the dynamical 
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Table I: /-sum rule tests for density and spin structure factors in 
BCS-BEC crossover, ujm is the maximum frequency. The expected 
value is 1.0 in units of hf€f/2. 
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structure factor for particle density can be written as the sum 
of five terms 

Sc — Sco + S MT c + S MT c + Sal + S co n, (6) 

where the first term on the right-hand-side denotes the "bare" 
term from Kq°. The remaining terms denote the corrections 
from the MT, and two AL diagrams and from the collective 
modes. 

In the spin channel, the spin-"magnetic field" interaction 
also contains effects associated with pairing fluctuations, but 
there is a significant difference. Here the AL diagrams are 
not included. This assertion can be verified by establishing 
numerical consistency with the /-sum rule (Q]). It is also, then, 
consistent with the equality in Eq. d2|. As a result, the full 
vertex is given by Ag = Xs + SA MT s + SA MT s . Here the bare 
spin-magnetic field interaction vertex is \s(P, P + Q) = 1. 
Diagrammatically these contributions are represented by the 
first two terms on the right hand side in Fig.Q] The full spin 
response function is given by xss(Q) = 

2^\ S (K,K~ Q)G{K)G{K - Q)A S (K, K - Q). (7) 

K 

We arrive at a final compact expression for the spin structure 
factor 

Ss = Sso + S MT s + S MJ s. (8) 

We now evaluate 5_(w,q). It can be shown that Sco = 
Sso- Moreover, the MT terms in the spin and density channel 
enter with reversed signs, as has been noted previously [10]: 
S MT c = — S MT s and 5 MT c = — 5 MT s. Together with the 
cancellation between MT and AL terms, we arrive at a very 
general form for 5_, defined in Eq.®. 

S- = S MT c + — Scon ■ (9) 

Importantly, S- (w, q) depends only on the order parameter; 
it vanishes above T c , as is trivially consistent with the sum rule 
in Eq. (0). For the strict BCS case, where the same behavior 
obtains, one can clearly see that this is not a consequence of 



Figure 3: (Color online) Dynamic density structure factor (solid 
lines) and spin structure factor (dashed lines) in BCS-BEC crossover 
for q — 0.5&;f. The black solid and dashed lines which show the be- 
havior slightly above T c , coincide with each other since Sc = Ss- 
The red lines indicate the behavior below T c where Sc 7^ Ss ■ 

the vanishing of the gap, but rather of true phase coherence. 
Above T c the same diagrams leading to the vanishing of the 
superfiuid density 111411 are responsible for the vanishing of 
S-. 

This discussion has been general; now we substitute our 
specific form for the self energy II 14fl to obtain 

<V q) or x ss ( W , q) = Y, [C 1 " /+ " /-) * 

p 

E+ + E+E_ - j.j + sgnA s 2 c - Ag g 
E+E_ w 2 -(E + +E_y [U } ~ ] 

E+ E_ E + E_ + - sgnA s 2 c + A p 2 g ] 

X E+E- u*-(E + -E_y I' (iU) 

Here E± — £p± q /2, f± — f(E±), and ui implicitly has a 
small imaginary part. The factor "sgn" is 1 and —1 in the 
density and spin channels respectively. Note that as a result 
of these sign changes in Eq. ( [Tol l, the dynamical correlations 
in the spin channel depend only on the square of the pair- 
ing gap A which involves the sum of the squares of A pg and 
A sc . This is very different from the particle density chan- 
nel which reflects the distinction between condensed (sc) and 
non-condensed (pg) pair contributions. 

We now establish that the sum rules (Eq.(fl]i) are satisfied 
along with Eq. (|2]). For illustration purposes, we choose a 
fixed momentum transfer q = Q.5}cf, and list the results of 
our sum rule checks for various scattering lengths Ufcl and 
two different temperatures (above and below T c ) in Table Q] 
Both sum rules are satisfied to within 10% or better. In the 
table ojm is the maximal frequency up to which we integrate 
in order to satisfy the sum rules. We found that both structure 
factors approach zero. On the BCS side, the structure factors 
decay extremely rapidly so that um need not be large, but 
on the BEC side of resonance, they exhibit a long and slowly 
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decaying tail. The characteristic size of this upper bound is 
of interest experimentally in measurements of the structure 
factor (which reflect the frequency integrated form of S(w, q)) 
such as in Ref. LL2J. We have also verified that the /-sum rule 
holds for large momentum transfers; at unitarity where our 
studies are most complete, the errors are smaller than 5%. 

Fig. |2] presents a plot of the the static structure factors 
5 c (q)/2 and S-(q) as a function the momentum transfer at 
unitarity and for two different temperatures below and above 
T c . The large q behavior at the lower T is consistent with the 
T = results in Ref. 9. At small q, all structure factors ap- 
proach zero as required by Eq.©. FigfJJ? is a central figure of 
this paper. Here we plot S- (q) as a function of temperature, 
which one can see behaves like an order parameter. The same 
result can be obtained for any fixed q and co and reflects the 
behavior implicit in Eq. (|9). That is, the difference between 
the density and spin channels arises only in the presence of 
a condensate. This behavior is also found in strict BCS the- 
ory as well as for any kpa. The vanishing of 5_ signals the 
vanishing of condensation. It should establish if an unknown 
Fermi gas is in the superfluid or normal phase. 

Fig- E] shows the particle density (solid lines) and spin re- 
sponse (dashed lines) in units of np/2ep as a function of 
uj/ep. We show temperatures both below and above T c for 
a fixed momentum transfer q = 0.5kp and for various kpa. 
Above T c , the magnetic and density response contributions 
coincide. Below T c , one sees from Fig. [3] that S c (u, q) has 
two main features, a collective mode peak, and a continuum 
of quasiparticle excitations which appears for u > 2 A, asso- 
ciated with the breaking of Cooper pairs induced by the per- 
turbation. The "collective" peak is broadened somewhat due 
to coupling to the continuum. As the system evolves from the 
BCS to the BEC regime (where the pairs are more difficult 
to break), this "collective" peak contains increasingly more of 
the total spectral weight contained in the /-sum rule. Above 
T c , the spectral weight in the former "collective" peak must be 
redistributed, leading to an increased contribution in the con- 
tinuum and to a new low frequency peak which is associated 
with thermally broken pairs. This peak appears at unitarity 
and on the BCS side of resonance for the considered range of 
temperatures. At a more formal level it arises from the con- 
tributions from the AL diagrams when /i > ^(f ) 2 above 
T c . These two continua nearly merge in the BCS regime in 
Fig- Ha); they are separated for the unitary and BEC cases. 

Turning now to the spin structure factor at low T (denoted 
by dashed lines in Fig. [3), it can be seen that there is no col- 
lective physics. The entire spectral weight of the /-sum rule is 
associated with the continuum, and this, in turn, only appears 
when u) > 2 A where the energy transfer is large enough to 
break pairs. Just as for the particle density response, we see 
that the low frequency spectral weight is highly suppressed. 
In the BEC regime, in Fig.[3|d), the spectral weight does not 
change significantly with T, since the temperatures consid- 
ered were not sufficient to break pairs. This behavior is con- 
sistent with that found earlier [11] on the BEC side of reso- 
nance. There is also a small low frequency peak, which gives 



a tiny contribution to the /-sum rule. 

In this paper we proposed a methodology for establishing 
phase coherence in neutral superfluids. This has been a long 
standing issue with one of the earliest such proposals based 
on noise correlations [17], although this has proved difficult 
to implement away from T = 0. Here, by contrast we pro- 
pose two photon Bragg experiments at T ^ which mea- 
sure the difference between the dynamical responses associ- 
ated with the particle and spin density, called 5_(w,q) (or 
fl flitl Sf±(ui, q)/2). The qualitative behavior shown in Fig- 
ure \2\b), also obtains for \S-(uj,q)\ for all q, ui. This sig- 
nature of phase coherent order, unlike others in the past based 
on Bragg scattering, requires looking at a single q,u),T. It 
should also apply to other crossover schemes which satisfy 
the conservation laws, to trapped gases and to optical lattices 
(at least) for frequencies below those of the band gaps. 
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